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The leading mean-field critical behaviour of </>|-theory is modified by multiplicative logarithmic corrections. We 
analyse these corrections both analytically and numerically. In particular we present a finite-size scaling theory 
for the Lee- Yang zeroes and temperature zeroes, both of which exhibit logarithmic corrections. On lattices from 
size 8 4 to 24 4 , Monte-Carlo cluster methods and multi-histogram techniques are used to determine the partition 
function zeroes closest to the critical point. Finite-size scaling behaviour is verified and the logarithmic corrections 
are found to be in good agreement with our analytical predictions. 
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1. INTRODUCTION 

The single component version of <fi theory in 
the d-dimensional Euclidean space-time contin- 
uum is defined by the Hamiltonian density 



n 



1 



;(V^) 2 



go 

4! 



+ %<f-H{x)<Kx) (1.1) 



where H(x) is the source for the fields <j){x). The 
lattice parameterization of the theory (in the ab- 
sence of a source field) is given by the action 
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Here the hopping parameter k and the quartic 
coefficient A correspond to the bare mass m and 
bare quartic coupling g respectively. The limit 
A — > 00 gives the Ising model. 

Above one dimension the discretized theory ex- 
hibits a phase transition (of second order) near 
which the continuum theory can be recovered. To 
remove the cutoff, it turns out that the quartic 
coupling has to be taken to the infra-red fixed 
point (IR FP) g R . The theory is believed to be 
trivial in d — 4 — although this has never been 
rigorously proved. This means that it is in the 
universality class of the theory of free bosonic 
fields. The leading (mean field) scaling behaviour 
is modified by logarithmic corrections, which are 
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linked to the triviality of the theory M. Their 
identification provides the primary motivation for 
this work. 

Logarithmic corrections to scaling in the infi- 
nite volume system have been studied in [|2J and 
||. Here we report on results for finite-size scal- 
ing (FSS) j|,[| of the <j) 4 theory which we have 
extended to four dimensions. Such finite size the- 
ories can be tested using non-perturbative (i.e., 
numerical) techniques. 

The usual statement of FSS is the following || : 
For any thermodynamic quantity Pl(k), mea- 
sured on a system of linear extent L and near 
criticality, 

Pl{k) .( L 
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where £00 (k) is the correlation length of the in- 
finite volume system. The usual justification 
for this formula is that L and £oo are the only 
length scales involved and hence their ratio, x = 
L/£oo(k), is the scaling variable. Until 1982 this 
statement had the status of a hyp othe sis. Then, 
Brczin Q4] succeeded in proving (L3) from the 
renormalization group (RG). An essential ingre- 
dient in this proof is that the running quartic cou- 
pling be approximated by its IR FP value g* R in 
the scaling region. Now, in d — 4 (in the per- 
turbative formulation at least), the IR FP of the 
Callan-Symanzik beta function is at the origin. 
The approximation above then leads to the mean 
field theory which predicts a phase transition even 



for a finit e system. For this reason FSS in the 
form Jl.g ) breaks down in d = 4. The intuitive 
justification given above is however a dimension 
independent argument. It is not clear, then, why 
it should fail in d = 4 while being valid for d < 4. 
We claim that the usual statement (|l.3| ) is, in 
fact, flawed and propose a modified FSS formula, 
valid in any dimension including four. 

2. THE PERTURB ATIVE RENORMA- 
LIZATION GROUP 

At some critical value, Wq c , of the bare mass, 
the renormalized theory is massless. Writing TOq 
as TOq c + 1, t becomes a measure of the deviation 
away from the massless theory. In the Ising ver- 
sion of the model, it is proportional to k — k c , k c 
being the critical hopping parameter. The gener- 
ating functional W[H, t] is defined by 
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C being a normalization constant. The function 
conjugate to H(x) is 
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If H is independent of x, (which we henceforth 
assume), then W is a function of its arguments. 
The generating functional T[M, t] of the one 
particle irreducible vertex functions is defined 
through the Legendre transformation 



T[M, t] + W[H, t}= j dxH(x)M(x), 



with 
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After isolating the divergences occurring in the 
Schwinger functions, one can write down the rela- 
tionship between the bare and renormalized the- 
ories. In order to be able to study the onset of 
criticality, in both the symmetric and the bro- 
ken phases, one first considers the massless renor- 
malized theory — renormalized at some arbitrary 
mass-scale parameter /i. Expanding in the re- 
duced mass t and in the conjugate function M, 



gives the renormalization group equation (RGE) 
for the massive theory in the critical region. Be- 
cause of the local nature of the renormalization 
group, the renormalization constants of the infi- 
nite volume theory render the finite volume the- 
ory finite too ||]. 

For a system of finite volume L d , with reduced 
temperature t and magnetization M, the above 
generating functional becomes the function 

r{t,M,g R ,n,L) 

in which g R represents the renormalized quartic 
coupling. The RGE expresses the invariance of 
the physics under a rescaling of fi. I.e., when 
the mass-scale fi is varied, t,M and g R respond 
in a way which is governed by the flow equations 
pLpfl . In four dimensions these flow equations can 
be solved perturbatively in g R . Rescaling fi to 
/z/L, and using dimensional analysis, gives the 
following solution of the RGE : 



T(t,M,g R ,l,L) 
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To determine how the running coupling on the 
right hand side of (2.8) couples to the remaining 
terms, perturbation theory must be applied to T 
itself. This gives |,§ 

r (t, M, g R , 1, L) = citM 2 (lnL)"5 

+c 2 M 4 (\n L)- 1 +c 3 t 2 (In L)i (2.9) 



where c\, . . . , C3 are constants. Applying (2.7) to 
this yields for the external field 



H(t,M,g R ,l,L)~ 

c 4 tM(lnL)-i + c 5 M 3 (In I)' 1 , 
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where, again, a and C5 are constants. 

These give for the free energy per unit volume 
in the presence of an external field 

W L (t,H) =c' 1 tM 2 (\nL)-^ 

+c' 2 M 4 (\nL)- 1 + c 3 t 2 (lnL)i, (2.11) 



c'j and c' 2 being consta nts an d M given by (|2.10 ) 
If H vanishes, then (2.1C) and (2.11) give 
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One could proceed directly from (2.11) or 
(2.12) to find the FSS formulae for thermody- 
namic observables. But it is more complete to 
study the partition function itself. This is entirely 
equivalent to the study of its zeroes. For fixed real 
t the zeroes in the complex h plane are called Lee- 
Yang zeroes M , and in the absence of an external 
field, the zeroes in t are called Fisher zeroes M. 
Their FSS properties below four dimensions was 
studied in pi. In this section, the corresponding 
FSS theory is presented for four dimensions where 
logarithmic corrections are manifest. 

The total free energy at the critical temper- 
ature in four dimensions in t he presence of an 
external field is given by ( 2.11 ) as 



L 4 (lnL)3#3. (2.13) 

The partition function is therefore 

Z L (t = 0, H) = Q (L 4 (In L)*H i) . (2.14) 



When the partition function is zero, solving for 
H gives 



Hj (xL- 3 (lnL)-- 
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where the constant of proportionality depends on 
the index j of the zero. This is the FSS formula 
for Lee- Yang zero es in four dimensions. 

If H vanishes, fl2.12 ) can be used in a similar 
way to show that the Fisher zeroes scale as 



tj cxL~ 2 (lnL) 
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Once the FSS behaviour of the partition func- 
tion zeroes has been found one can easily find 
the corresponding behaviour for thermodynamic 
functions by expressing them in terms of the ze- 
roes. These considerations give for the zero field 
magnetic susceptibility and specific heat 
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3. NON-PERTURBATIVE ANALYSIS 
OF FINITE SIZE SCALING 

The Swendsen-Wang cluster algorithm jlQ] was 
applied to the Ising version of the theory on lat- 
tices of sizes 8 4 to 24 4 . 

In an external field h (= kH), the partition 
function can be written as 



Z(k, h) 



where 
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and the spectral density p(S, M) is the relative 
weight of configurations having given values of S 
and M. The 'multihistogram' method |11|] was 
used to combine histograms determined at vari- 
ous values of k. This provides an optimal esti- 
mator for the spectral density and allows one to 
construct Z(k, h) in the complex neighbourhood 
of the critical point. A Newton-Raphson algo- 
rithm was used to determine nearby zeroes. 

The leading (power law) FSS behaviour of the 
Lee- Yang and Fisher zeroes was found to be 
slightly deviant from the mean field predictions. 
These deviations find their explanation in the 
presence of logarithmic corrections. To isolate 
these corrections, in the case of Fisher zeroes, we 
plot in fig. la ln(L 2 ImKi) versus ln(lnL). The 
negative slope is in good agreement with the scal- 
ing prediction of — |. In fact, a fit to all five 
points gives a slope —0.217(12). Excluding the 
point corresponding to L — 8 gives a slope of 
-0.21(4). The solid line is the best fit to the re- 
maining point s assu ming the theoretical predic- 
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We may now determine k c from \kj — n c 
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Using the first Fisher zeroes, we 



find k c ~ 0.149703(15) in good agreement with 
the value 0.149668(30) from high temperature ex- 
pansions jjjj. 

To identify the logarithmic corrections for the 
Lee- Yang zeroes, we plot in fig. lb ln(L 3 Im/ii) 
against In (In L). A best fit to all five points gives 
a slope of —0.204(9) which compares well with 
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the theoretical prediction of — \ from ( 2.15 ). 
eluding the smallest lattice, a fit to the remaining 
four points gives a slope —0.22(3). The solid line 
in fig. lb is the best fit to the last four points with 
given slope —j. 



Figure 1. Logarithmic corrections to FSS of (a) 
Fisher zeroes and (b) Lee- Yang zeroes. 



4. CONCLUSIONS 

A finite size scaling theory has been developed 
for the single component </> 4 theory in d = 4 
dimensions. Emphasis has been placed on log- 
arithmic corrections to the mean field predic- 
tions. This has been checked non-perturbatively 
using high precision numerical methods, and good 
agreement is found. 

FSS formulae for other thermodynamic func- 
tions are also given. These exhibit logarithmic 
corrections too. The FSS formula for the correla- 
tion length of a four dimensional system also in- 
volves logarithmic corrections. This was derived 



by Brezin |4J| for a system of extent L in all direc- 
tions. At the infinite volume critical point k — k c , 



£l(k c ) oc L(lnL)" 
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This suggests that the FSS variable should be 
&(k c ) L(\nL)i 



£oo(k) 
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in four dimensions |l3| . Indeed, replacing t he s cal- 
ing variable, x. of the right hand side of fll.3| ) by 
the ratio £l(k c )/£oo(k) is sufficient to recover all 
the FSS formulae presented here while still being 
correct in d < 4 dimensions. We suggest that this 
modified FSS hypothesis is the more appropriate 
one. 
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